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In this work, we present a systematic derivation of the distribution of eigenfrequencies for oscil- 
lations of the ground state of a repulsive Bose-Einstein condensate in the semi-classical (Thomas- 
Fermi) limit. Our calculations are performed in 1-, 2- and 3-dimensional settings. Connections with 
the earlier work of Stringari, with numerical computations, and with theoretical expectations for 
invariant frequencies based on symmetry principles are also given. 

I. INTRODUCTION 

Bose-Einstein condensation (BEC) is one of the most exciting achievements within the physics community 
in the last two decades. Its experimental realization in 1995, by two experimental groups using vapors of Rb 
[l[ and Na Q marked the formation of a new state of matter consisting of a cloud of atoms within the same 
quantum state, creating a giant matter wave. However, in addition to its impact on the physical side, this 
development had a significant influence on mathematical studies of such Bose-Einstein condensates (BECs) 
@j EL HI, @, 0]- I n considering typical BEC experiments and in exploring the unprecedented control of the 
condensates through magnetic and optical "knobs" , a mean-field theory is applied to reduce the quantum 
many-atom description to a scalar nonlinear Gross-Pitaevskii equation (GPE). This is a variant of the famous 
nonlinear Schrodinger (NLS) equation [H, Q of the form: 

d^> ti 2 

ih— = -—X7 2 y + g\y\H + V cxt (r)>S>, (1) 

where = ^(r, t) is the BEC wavefunction (the atomic density is proportional to |^(r, t)\ 2 ), V 2 is the Laplacian 
in r = (.x, y, z), m is the atomic mass, the prefactor g is proportional to the atomic scattering length (e.g. g > 
for Rb and Na, while g < for Li atoms), and Kxt(r) is the external potential for magnetic or optical traps. 

The NLS equation is a well-established model in applications in optical and plasma physics as well as in fluid 
mechanics, where it emerges out of entirely different physical considerations [81, 95]. In particular, for instance, 
in optics, it emerges due to the so-called Kerr effect, whereby the material refractive index depends linearly on 
the intensity of incident light. The widespread use of the NLS equation stems from the fact that it describes, 
to the lowest order, the nonlinear dynamics of envelope waves. 

One of the particularly desirable features of GPE is that the external potential V ex t(r) can assume a multi- 
plicity of forms, based on the type of trapping used to confine the atoms. Arguably, however, the most typical 
magnetic trapping imposes a parabolic potential [5j, [lCj 



V^^^iu^+uy+^z 2 ), (2) 



where, in general, the trap frequencies u x ,y,z are different. In what follows, however, for simplicity, we will 
restrict our consideration to the isotropic case of equal frequencies along the different 

directions. 

If the wave function is decomposed according to 

tf(r,t) =e~ ltlt/h U(r), 

then U(r) solves the stationary GPE with the chemical potential /i. One of the most extensively discussed 
limits in the case of self-repulsive nonlinearity g > and in the presence of the parabolic potential ([2]) , is the 
limit fx — ► oo. This limit is referred to as the Thomas-Fermi limit. 
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If the kinetic (Laplacian) term is neglected, the approximate ground state solution is obtained in the form 

U(r) = (maxLg-^-KxtMXO]) 172 . 

For this limit, the seminal work of Stringari (ll| suggested a computation of the corresponding eigenfrequencies 
of oscillations of perturbations around the ground state of the system, using a hydrodynamic approach. This 
approach has become popular in the physics literature for more complicated problems involving anisotropic 
traps and dipole-dipole interactions [T3 |. 

The aim of the present work is to derive these eigenfrequencies systematically not only in the 3-dimensional 
context, but also in the 2-dimensional and 1-dimensional cases. We will relate these eigenfrequencies to the 
eigenvalues discussed in the recent work [14] in the 1-dimensional setting. The relevant eigenfrequencies of 
the perturbations around the ground state are then directly compared with numerical computations in the 
2-dimcnsional case. 

In the numerical computations, the eigenfrequencies are obtained systematically as a function of the chemical 
potential \x starting from the low-amplitude limit (when the ground state is approximately that of the parabolic 
potential) all the way to the large-chemical potential. Earlier, these eigenfrequencies were approximated 
numerically near the low-amplitude limit by Zezyulin et al. fl5j in the 1-dimensional case and by Zezyulin [l6j 
in the 2-dimensional case. Our numerical computations also allow us to identify eigenfrequencies that remain 
invariant under changes in \x and to connect them to underlying symmetries of the GPE. 

Our presentation will be structured as follows. In section 2, we present the mathematical setup of the problem. 
In section 3, we compute its corresponding linearization eigenvalues (around the Thomas-Fermi ground state) 
in the 1-, 2- and 3-dimensional settings. In section 4, we compare these results to direct numerical computations 
in the 2-dimensional case. Note that the 1-dimensional case was considered in some detail in our earlier work 
[l?]]. Lastly, a brief summary of our findings and some interesting directions for future study are offered in 
section 5. 



II. MATHEMATICAL SETUP 



Using rescaling of variables, one can normalize the GPE (JTJ into two equivalent forms. One form corresponds 
to the semi-classical limit and it arises if h = e, m = i, g = 1, uj = 2, and /i = 1, or equivalently, in the form 

ieut + e 2 V 2 u + (1 - \x\ 2 - \u\ 2 )u = 0, (3) 

where u(x,t) : R d x M — * C is a wave function, V 2 = d%. + ... + d 2 d is the Laplacian operator in d spatial 
dimensions, and e is a small parameter. On the other hand, if 

v = fi 1/2 u, £ = (2fi) 1/2 x, t = 2t, (4) 

then equation ([3]) can be translated to the form 

iv T =-^lv+^\ 2 V+\v\ 2 V~ fiV, (5) 

that corresponds to the GPE (JTJ with h=l,m = l, g = l,u=l, and [i = The semi-classical limit e — » 
corresponds to the Thomas-Fermi limit fi — ► oo. 

Let % be a real positive solution of the stationary problem 

e 2 V 2 r; £ + (1 - |a;| 2 - r) 2 )r)e = 0, x e M d . (6) 

According to Gallo & Pelinovsky [l8| , for any sufficiently small e > there exists a smooth radially symmetric 
solution r\ e 6 C 00 (]R d ) that decays to zero as |x| — > oo faster than any exponential function. This solution 
converges pointwise as e — > to the compact Thomas-Fermi cloud 

,. / (1-M 2 ) 1/2 , for \x\ < 1, (r . 

?7 := urn r? £ = ^ V r , (7) 

' e^o ' \ 0, for \x\ > 1. w 
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The solution r\ e with the properties above is generally referred to as the ground state of the Gross-Pitaevskii 
equation ([3|). 

The spectral stability problem (often referred to as the Bogolyubov-de Gennes problem in the context of 
BECs) for the ground state n £ is written as the eigenvalue problem 

L + u = —Xew, L_w = Xeu, (8) 

associated with the two Schrodinger operators 

J L+ = -£ 2 V 2 + |x| 2 +3r/ 2 -l, 
\ L_ = -e 2 V 2 + jirj 2 + n 2 - 1. 

A naive approximation of the eigenvalues A of the spectral stability problem ([8]) arises if we replace rj 2 by 
rj 2 . Because L + is invertible, the eigenvalue problem can then be written in the form 

(-e 2 V 2 + |x| 2 + 7 7 2 -l)u; = 7£ 2 (-e 2 V 2 + |2;| 2 + 3r; 2 -l) _1 i i ;, x E R d , (9) 

where 7 = —A 2 . The formal limit e — > gives a restricted problem in the unit ball 

LI: -2(1 - \x\ 2 )V 2 w = 710, x e B = {x e R d : \x\ < 1}, (10) 

subject to the Dirichlct boundary condition on the sphere |x| = 1. Convergence of eigenvalues of to 
eigenvalues of the limiting problem (fT0| was rigorously justified by Gallo & Pelinovsky [14| in one spatial 
dimension d = 1. 

Because L + is invertible for any small e > 0, the original eigenvalue problem ((5|) can also be written in the 
form 

(-e 2 V 2 + \x\ 2 +r)l - 1) w = 7 e 2 (-e 2 V 2 + \x\ 2 +3r/ 2 - l) _1 w, xeR d . (11) 

Because 

\x\ 2 + r, 2 -l = ^, 

Ve 

the formal limit e — > gives now a different problem in the unit ball 

LII: -2(1 - Id 2 ) ( V 2 w - ^-^-w] = jw, xeB , (12) 

V m J 

subject to Dirichlet boundary conditions on the sphere \x\ — 1. Justification of convergence of eigenvalues of 
(jlip to eigenvalues of the limiting problem \12\ is still an open problem in analysis. 
The limiting eigenvalue problem (|10[) can be written in the vector form 

LI: -2V(1- |x| 2 )Vv = 7 v, x € B , (13) 

where v = Vw G M. d . On the other hand, the limiting eigenvalue problem (|12[) can be rewritten in the equivalent 
scalar form 

LII: -2V(1 - \x\ 2 )Vv = jv, x e B , (14) 

where v = and t)q is given by ([7]). It was exactly the representation (TT4"|) of the limiting eigenvalue problem 
LII, which was derived by Stringari [ll[ from the hydrodynamical formulation of the Gross-Pitaevskii equation 
§3§ in three dimensions d = 3. 

Comparison of the two representations (| 13() and (|14p implies that the two limiting eigenvalue problems have 
identical nonzero eigenvalues in the space of one dimension d = 1 but may have different nonzero eigenvalues 
for d 2. We will show in the next section that it is exactly the case. We will illustrate numerically for d = 2 
that the eigenvalues of the second limiting problem (TT4"]) are detected in the limit e — ► from the eigenvalues 
of the original problem ([8]). 
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III. EIGENVALUES OF THE LIMITING PROBLEMS 



Case d = 1: Both representations (fT5|) and (fT4")) of the limiting eigenvalue problems LI and LI I reduce to 
the Legendre equation 

(x 2 - l)v"(x) + 2xv'(x) = ]pv{x), (15) 

For LI given by (fT0|) . the correspondence of eigenfunctions is v{x) = w'{x). The only nonsingular solutions 
of this equation at the regular singular points x = ±1 are Legendre polynomials v(x) £ {P n (x)} n ^o, which 
correspond to eigenvalues 7 £ {2n(n + l)}„^o- The zero eigenvalue must be excluded from the set since it 
corresponds to w(x) = x, which violates Dirichlet boundary conditions at x = ±1 for w(x). On the other hand, 
all nonzero eigenvalues are present because the corresponding eigenfunction w(x) = C n+1 (x) constructed from 
v(x) = P n (x) thanks to identities 8.936, 8.938, and 8.939 in [ijf 

±C-l{ 2 (x)=CV 2 (x) = P n (x), 
also satisfies the Dirichlet boundary conditions thanks to the identity 

For LI I given by fl2")) . the correspondence of eigenfunctions is v(x) = 1^ a ■ The zero eigenvalue should 

now be included for n = 0, since the eigenfunction w(x) = \J\ — x 2 corresponds to the ground state r\ e in the 
limit £ — > 0, which is known to be the eigenfunction of operator L_. All nonzero eigenvalues are the same as 
for the limiting problem (|10p but the eigenfunctions are now different. For the eigenvalue 7 = 2n(n + 1), the 
eigenfunction is w(x) = yl — x 2 P n (x). It should be noted here that the obtained eigenvalue distribution was 
numerically examined in [ItJ and was found to be in very good agreement with the true eigenvalues of system 
©. 

Case d = 2: We consider the limiting eigenvalue problem LI I in the form (fT4"]l and use the polar coordinates 

x = r cos(0), . _ n r „ _ , 
• (J r > 0, € 0,2tt . 
y = rsm(0), ' L ' J 

After the separation of variables v(r, 0) = V (r)e l for m £ Z, we obtain an infinite set of eigenvalue problems 
for amplitudes of cylindrical harmonics 



(1 - r 2 ) (v"(r) + -V'(r) - ^-V(r)) + 2rV'(r) - \~{V(r), r £ (0, 1). 
\ r r z J 2 



(16) 



Let m ^ 0. We are looking for solutions of equation (|16p which behave like V(f) ~ r m as r — > 0. Let us 
transform (fl~6)) to a hypergeometric equation with the substitution V(r) = r m F(z), z = r 2 . Direct computations 
show that F(z) solves 



z(l - z)F"(z) + (1 + m - (2 + m)z)F'(z) + Q7 - im^j 



F(z) = 0, 2 6(0,1). 



A nonsingular solution at z = is the hypergeometric function F{z) = T{a, b, c; z) where 

vn 7 

c = 1 + m, a + 6 = 1 + m, a6= . 

' ' 2 8 

Because a + b — c = 0, the hypergeometric function is singular at z = 1 unless it becomes a polynomial for 
a = — fc with an integer fc 0. The eigenvalues of the limiting problem (I14p are then given by 

7 e {tSWo,^, 7Si = 4(m + 2fe(l + m) + 2fc 2 ). (17) 
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Let us now consider the limiting eigenvalue problem LI in the form (jlOp and use the same polar coordinates. 
The corresponding eigenvalue problem is 

^- 2 >(^ + ^4S0=™ <,e[0,24 (18, 

Let w(r, 0) — V(r)e lmS for m £ Z and obtain an infinite set of eigenvalue problems for cylindrical harmonics 
-(l-r 2 )(v'\r) + ~V'(r)-^V(r)^ =~jV(r), re (0,1). (19) 

Let m ^ 0. Equation (JT9J) can also be transformed to a hypergeometric equation after the substitution 
V(r) — r m F(z), z — r 2 . Direct computations show that F(z) solves 

z(l - z)F"(z) + (1 +m)(l - z)F'(z) + -~fF{z) = 0, z e (0, 1). 

o 

A nonsingular solution at z = is the hypergeometric function F(z) = T{a, b, c; z) where 

7 

c = 1 + m, a + 6 = to, ab = — — . 

8 

Because a + b — c = — 1, the hypergeometric function is bounded at z = 1. However, we need = to 

satisfy the Dirichlet boundary conditions for w(r, 6) at r = 1. From this condition, -F(z) has to be a polynomial 
(or F'(z) and F"(z) are singular at z = 1 with .F(l) 7^ 0). The polynomial arises for a = —k with an integer 
k 0. The eigenvalues of the limiting problem (|10[) are then given by 

T6{t5Wo^o, 7i 2) fc = 8fc(m + fc). (20) 

Comparison of (fT7|) and (j2"0)) show that 7^ = 7^ for all k ^ 0, but the sets {T^'jm^oi^o and 
{7m,fc}™^o,fe>o are different. For instance, 7}^ = 4 is not present in the set {j^ k } m ^oM^o- 

Case d = 3: We consider the limiting eigenvalue problem LI I in the form (fT4|) and use the spherical 
coordinates 

x = r cos(0) cos((/p), 

y = rsin(0)cos(<p), r ^ 0, 0e[O,27r], ^g[0,tt]. 
z = r sin(<p), 

After the separation of variables v(r, 9, cp) = V(r)Yi >m (6, <p) for m € Z and / G N, where Y^ m {6, <p) are spherical 
harmonics, we obtain an infinite set of eigenvalue problems for amplitudes of the spherical harmonics: 

-(1-r 2 ) (V"(r) + -V'{r)- l -^^-V(r)] + 2rV'(r) = ljV(r), re (0,1). (21) 
\ r r z J 2 

Using a similar reduction V(r) — r l F(z), z = r 2 to the hypergeometric equation, we obtain the eigenvalues of 
the limiting problem (|14[) in the form 

7 G {-Y$h>o,k>o, 7$ = 4(Z + 3fc + 2fcZ + 2fc 2 ). (22) 
This distribution was obtained by Stringari [llT | from the balance of the leading powers in polynomial solutions 

ofinsj. 

Using the same algorithm, the eigenvalues of the limiting problem LI are found in the form 

7 G {7$}j>o,k>o, 7$ - 4fc(l + 21 + 2k). (23) 
This distribution is different from (|22[) . In particular, it does not include eigenvalue 7^ = 4. 
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IV. NUMERICAL RESULTS 



As indicated above in the one-dimensional case, good agreement was observed between the predicted Thomas 



Fermi limit spectrum and the numerical computations of 17]; for this reason, we now turn our attenion to the 
two-dimensional case. Eigenvalues of the original spectral problem ([8]) for d = 2 are computed numerically 
and shown on Figure [TJ (solid lines) together with the limiting eigenvalues (jTT)) of the reduced spectral problem 
LI I (dash-dotted lines). Notice that the results are presented in the context of the rescaled variant of the 
Gross-Pitaevskii equation ^ commonly used in the physical literature, illustrating the relevant eigenvalues as 
a function of the chemical potential fj, = 1/ (2s) . 

The ground state r\ e exists for any s < \ (i.e., /i > 1) and the limit e — + h can be obtained via small- 
amplitude bifurcation theory [lj|. All eigenvalues 7 = — A 2 in the spectral problem {8]| in this limit occur at 
the integers A(n + m) 2 with n, m ^ and the multiplicity of the eigenvalue 7 = 4(n + m) 2 is n + m + 1. When 
e < i, this degeneracy is broken and all eigenvalues become smaller as e gets smaller (or /1 increases) besides 
the double eigenvalue 7 = 4 and the simple eigenvalue 7 = 16. Notice that the eigenvalues 7 are related to the 
eigenfrequencies lo on Figure [TJ by w — (accounting for the time rescaling r = 2t). 

Persistence of e-independent eigenvalues 7 = 4 and 7 = 16 of the spectral problem ([8]) is explained by the 
symmetries of the Gross-Pitaevskii equation One symmetry is given by the explicit transformation of 
solutions 

u(x 7 y,t) = e^ x+ls ^y + ^u(x,y,t), x = x - q(t), y = y-k(t). (24) 

If u(x , y, t) is a solution of equation ([3]) rewritten in tilded variables and (p, s, q, fc, lo) satisfy 

q = 2ep, ep + 2q = 0, 1 
k = 2es, ss + 2k = 0, - = -2^ + fcs ), 

then u(x,y,t) is also a solution of equation ([3]). Therefore, both q and k satisfy the linear oscillator equations 
with eigenvalue 7 = 4, which gives the double degeneracy of eigenfrequency lo = 1 in Figure [TJ 

The other symmetry of the Gross-Pitaevskii equation ([3]) with d = 2 is given by the conformal transformation 

u(x,y,t) =a(t)u(x,y,t)e lc(t)x2+lc{t)v2+luj{t \ x = a(t)x, y = a(t)y, t = b(t), (25) 

where (a, b, c, lo) satisfy the first-order differential equations 

b = a 2 , euj = l-a 2 , a + Aeac = 0, ec + Ac 2 e 2 + 1 = a 4 . 

Excluding c and denoting a(t) = z _1 (i), we obtained the nonlinear oscillator equation for z{t): 

z + Az{l - z~ 4 ) = 0. 

There is a unique critical point z = 1 and it is a center with eigenvalue 7 = 16, corresponding to the eigenfre- 
quency lo = 2 in Figure [TJ 

The two symmetries (|24|) and (|25|) explain the e-independent eigenfrequencies lo on Figure [TJ On the other 
hand, the figure shows that all eigenvalues 7 approach to the limiting eigenvalues (|17p as e — > (i.e., as /i — ► 00). 
This output confirms the robustness of the asymptotic distributions presented herein. 

It is worth noting that in all the cases shown on Figure [TJ the eigenvalues have been confirmed including also 
their multiplicities. For instance the eigenfrequency associated with lo = 2 (7 = 16) is associated with two 
eigenvalues in the set (fTT|) . namely with k = 1 and m = 0, as well as with k = and m = 4. One of these 
corresponds to the conformal symmetry (|2"5|) , while the other one can be observed on Figure [TJ to asymptote 
to lo = 2 in the large-// limit, as expected. 



V. CONCLUSIONS AND FUTURE DIRECTIONS 



In the present work, we offered a systematic approach towards identifying the eigenfrequencies of oscilla- 
tions of the perturbations around the ground state of a Bose-Einstein condensate in an arbitrary number of 
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FIG. 1: Eigenfrequencies of the two-dimensional spectral problem |(5J (solid blue lines) and the limiting eigenvalues 
of the reduced problem (|12[) (dash-dotted red lines). Notice that for practical numerical reasons the computation was 
performed for a parabolic trap of frequency fl = 0.3 and the corresponding eigenfrequencies and chemical potentials 
were appropriately rescaled by fi. 

dimensions (our calculations were given in 1-, 2- and 3-dimensions) . This spectrum is important because it 
corresponds to the excitations that can be (and have been) experimentally observed once the condensate is 
perturbed appropriately. 

Part of the rationale for attempting to understand the details of this spectrum is that when fundamental 
nonlinear excitations are additionally considered on top of the ground state, then the spectrum contains both 
a "ghost" of the spectrum of the ground state and the so-called negative energy modes that pertain to the 
nonlinear excitation itself. Relevant examples of this sort can be found both for the case of one-dimensional 
dark solitons (and multi-solitons) as analyzed in (20| and in the case of two-dimensional vortices, as examined 
in plj . It is then of particular interest to try to understand eigenfrequencies of excitations of these structures in 
the Thomas-Fermi limit, as well as those of their three-dimensional generalizations bearing line- or ring- vortices. 
Such studies would be especially interesting for future works. 

Acknowledgments: PGK is partially supported by NSF-DMS-0349023 (CAREER), NSF-DMS-0806762 
and the Alexander-von-Humboldt Foundation. DEP is supported by the NSERC grant. 



[1] M.H.J. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, "Observation of Bose-Einstein condensation in a 
dilute atomic vapor", Science 269, 198-201 (1995). 

[2] K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Durfee, D.M. Kurn and W. Ketterle, "Bose- 
Einstein condensation in a gas of sodium atoms", Phys. Rev. Lett. 75, 3969-3973 (1995). 

[3] C.J. Pethick and H. Smith, Bose-Einstein condensation in dilute gases, Cambridge University Press (Cambridge, 
2002). 

[4] L.P. Pitaevskii and S. Stringari, Bose-Einstein Condensation, Oxford University Press (Oxford, 2003). 

[5] F. Dalfovo, S. Giorgini, L.P. Pitaevskii and S. Stringari, "Theory of Bose-Einstein condensation in trapped gases", 
Rev. Mod. Phys. 71, 463-512 (1999). 

[6] P.G. Kevrekidis, D.J. Frantzeskakis, and R. Carretero-Gonzalez (eds.). Emergent Nonlinear Phenomena in Bose- 
Einstein Condensates: Theory and Experiment. Springer Series on Atomic, Optical, and Plasma Physics 45 



8 



(Springer, Heidelberg, 2008). 

[7] R. Carretero-Gonzalez, D.J. Frantzeskakis, and P.G. Kevrekidis. "Nonlinear Waves in Bose-Einstein Condensates: 

Physical Relevance and Mathematical Techniques", Nonlinearity 21, R139-R202 (2008). 
[8] C. Sulem and P.L. Sulem, The Nonlinear Schrodinger Equation, Springer- Verlag (New York, 1999). 
[9] M.J. Ablowitz, B. Prinari and A.D. Trubatch, Discrete and Continuous Nonlinear Schrodinger Systems, Cambridge 
University Press (Cambridge, 2004). 

[10] P.G. Kevrekidis and D.J. Frantzeskakis, "Pattern forming dynamical instabilities of Bose-Einstein condensates", 
Mod Phys. Lett. B 18, 173-202 (2004). 

[11] S. Stringari, "Collective excitations of a trapped Bose-condensed gas", Phys. Rev. Lett. 77, 2360-2363 (1996) 

[12] M. Fliesser, A. Csordas, P. Szepfalusy, and R. Graham, "Hydrodynamic excitations of Bose condensates in 
anisotropic traps", Phys. Rev. A 56, R2533-R2536 (1997) 

[13] C. Eberlein, S. Giovanazzi, and D.H.J. O'Dell, "Exact solution of the Thomas-Fermi equation for a trapped Bose- 
Einstein condensate with dipole-dipole interactions", Phys. Rev. A 71, 033618 (2005) 

[14] C. Gallo and D. Pelinovsky, "Eigenvalues of a nonlinear ground state in the Thomas-Fermi approximation", J. 
Math. Anal. Appl. 355, 495-526 (2009) 

[15] DA. Zezyulin, G.L. Alfimov, V.V. Konotop, and V.M. Perez-Garci'a, "Stability of excited states of a Bose-Einstein 
condensate in an anharmonic trap", Phys. Rev. A 78, 013606 (2008) 

[16] DA. Zezyulin, "Stability of two-dimensional radial excited states of a Bose-Einstein condensate in an anharmonic 
trap", Phys. Rev. A 79, 033622 (2009) 

[17] D.E. Pelinovsky and P.G. Kevrekidis, "Periodic oscillations of dark solitons in parabolic potentials", Cont. Math. 
473, 159-179 (2008) 

[18] C. Gallo and D. Pelinovsky, "On the Thomas-Fermi ground state in a radially symmetric parabolic trap", 
larXiv:0911.3"9l3l f2009) 

[19] I.S. Gradshteyn and I.M. Ryzhik, Table of integrals, series and products, 6th edition, (Academic Press, 2005) 
[20] G. Theocharis, A. Weller, J. P. Ronzheimer, C. Gross, M.K. Oberthaler, P.G. Kevrekidis and D.J. Frantzeskakis, 

"Multiple atomic dark solitons in cigar-shaped Bose-Einstein condensates", arXiv:0909.2122 
[21] S. Middelkamp, P.G. Kevrekidis, D.J. Frantzeskakis, R. Carretero-Gonzalez and P. Schmelcher, "Anomalous modes 

and matter-wave vortices in the presence of collisional inhomogeneities and finite temperature", arXiv:0911.3308 



